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• We present the new preliminary MRST parton distributions at NLO and NNLO. 

\J^) ' The analysis includes some new data and there is an improvement in the theoretical 

f**^ treatment at NNLO. Essentially complete NNLO partons are presented for the first 

, time, together with uncertainties. 

There are a number of reasons for an update of the MRST parton 
distributions. Firstly, there are new data to be included: NuTeV data 1 
on F£*(x, Q 2 ) and F^(x, Q 2 ) replacing CCFR 2 ; new CDFII high-^ T jet 
I data 3 (only compared, not yet fit); and we now include direct high- a: data 

on Fl(x,Q 2 ). There are also major changes in the theory: an implemen- 
tation of a new heavy flavour VFNS 4 , particularly at NNLO; and the in- 
clusion of NNLO corrections 5 to the Drell-Yan cross-sections. This leads 
to some important changes as NLO — > NNLO. The most important change 
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compared to the previous NNLO partons 6 , which already used the exact 
splitting functions 7 , is the new VFNS which leads to a significant change 
in the gluon and heavy quarks. Moreover, due to the NNLO procedure 
being essentially complete we now examine the uncertainties on the NNLO 
partons. In general the size of the uncertainties due to experimental errors 
is similar to that at NLO 8 . There is more work to do in order to estimate 
the theoretical uncertainty, which is certainly important in some regions 9 . 

We first consider the new data in the fit. The NuTeV structure func- 
tion data are not completely compatible with the older CCFR data. The 
main source of the discrepancy is in the calibration of the magnetic field 
map of the muon spectrometer, i.e. in the muon energy scale. However, 
the previous parton distribution fits were perfectly compatible with the 
CCFR data using an EMC inspired Q 2 -independent nuclear correction 10 
R. This correction is far too large for the new NuTeV data. The high- a; 
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region is completely dominated by the valence quarks for both F^'" (x, Q 2 ) 
and F^' v (x,Q 2 ). These are well known from fixed target F^' p (a;, Q 2 ) and 
F fl,d (x,Q 2 ). In order to fit the NuTeV data we try a reduced correction 
factor R eff = 1 + A* (R - 1). The best fit is for A = 0.2 and the previous 
nuclear correction is clearly ruled out. Hence, the NuTeV data imply a 
nuclear correction which is different for neutrinos than for charged leptons. 
However, recent CHORUS 11 data are in much better agreement with the 
CCFR data than the NuTeV data. Also, the partons in the region of high 
nuclear correction are already well determined. It may be appropriate to 
cut the nuclear target data in this region. The important information that 
neutrino DIS gives on the flavour composition of the proton is in the region 
x < 0.3, where the nuclear corrections are not so large or uncertain. 
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Figure 1. Comparison of the NLO up distribution with the NNLO up distribution, 
concentrating on small x (left) and high x (right). 

The comparison of the most recent MRST partons 6 can be seen along- 
side the Tevatron jet data 3 . The prediction clearly matches quite well, and 
it is clear that the new data still require the large high- a; gluon in the MS 
scheme. The fit to the fixed target data on F L (x, Q 2 ) also prefers the larger 
gluon since the data are generally larger than NLO or NNLO 12 , and a large 
coupling (and/or higher twist contributions) is needed. 

The change in the up distribution when going from NLO to NNLO is 
shown in Fig. 1. At small x the effect of the coefficient functions, partic- 
ularly C2, g (x,Q 2 ), is important and the difference between the NLO and 
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the NNLO distribution is greater than the uncertainty in each calculated 
using the Hessian approach 13 . At large x the coefficient functions are again 
important - C 2 q (x) ~ (ln 3 (l — x)/(l — x)) + and the difference between 
NLO and NNLO is again larger than the uncertainty in each. There is no 
real change from the MRST2004NNLO partons for the light quarks. At 
small x the effect of the splitting functions is important, particularly from 
P 2 g (x,Q 2 ), which has a positive ln(l/x)/x contribution. This affects the 
gluon distribution via the fit to dF 2 {x,Q 2 )/d\iiQ 2 , and the NNLO gluon 
is smaller at very low x than the NLO gluon. 
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Figure 2. Comparison of the NNLO MRST gluon distribution with the previous ap- 
proximate NNLO distribution (left) and the NNLO distribution of Alekhin (right). 

At NNLO heavy flavour no longer evolves from zero at /i 2 = m 2 , i.e. 
(c+c)(x,ml) — A 2 Hg {m 2 c )®g{m 2 ). In practice it starts from a negative value 
since the matrix element is negative at small x. The increased evolution 
from the NNLO splitting function allows the charm distribution to catch 
up partially with respect to that at NLO, which starts from zero at m 2 but 
it always lags a little at higher Q 2 . The correct NNLO charm is smaller 
than the approximate MRST2004 distribution which turned on from zero. 
This correction in the charm procedure also affects the gluon compared to 
the MRST2004 NNLO partons, Fig. 2, and the change is greater than the 
uncertainty in some places. The correct heavy flavour treatment is vital. 

At NNLO the Drell-Yan corrections 5 are significant. There is an en- 
hancement at high xf — x\—X2 due to large logarithms, which is similar to 
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the ln(l — x) enhancement in structure functions. The NLO correction is 
large and the NNLO corrections are 10% or more. The quality of the fit to 
E866 Drell-Yan production 14 in proton-proton collisions is x 2 = 223/174 
at NLO and x 2 = 240/174 at NNLO. The scatter of points is large and 
a x 2 ~ 220 is the best possible. The quality of the fit is good, as seen in 
Fig. 3. It is worse for proton-deuteron data. The positive correction at 
NNLO requires the data normalization to be 110% (103% at NLO), there 
being little freedom since the sea quarks for x < 0.1 and the valence quarks 
are already well determined by structure function data. The normalization 
uncertainty on the data is 6.5%, and a change of 10% is a little surprising. 

The quality of the full fit at NLO is X 2 = 2406/2287 and at NNLO is 
X 2 = 2366/2287. NNLO is fairly consistently better than NLO. There is 
a definite tendency for as(M§) to increase with all changes, both the new 
data and the improved theoretical treatment. At NLO as(M|) = 0.121 and 
at NNLO as(Mz) = 0.119. Although the fit is generally good, particularly 
at NNLO, there is some room for improvement, and the data would prefer 
a little more gluon at both high and moderate x. 
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Figure 3. Comparison of the NLO MRST Drell-Yan cross-sections with the data. 

We compare with the only other NNLO partons available, those of 
Alckhin 15 (there is nothing at present from CTEQ). We have a much larger 
as(Mz), i-e. as(M§) = 0.119 compared to 0.114. There is not much dif- 
ference in high- a; valence quarks, except that explained by the difference 
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in as(Mz). There are differences in the low-x sea quarks but these are 
dominated by differences in flavour treatments of u — d and s(x,Q 2 ). The 
gluon distribution difference at small x is seen in Fig. 3, and is much bigger 
than the uncertainties. This is due to the heavy flavour treatments, which 
we have already shown to be important as well as to differences in the data 
fit and in as(M§ ). The gluons also differ a great deal at high x, where they 
are determined by the Tevatron jet data 16 for MRST, the comparison now 
being excellent 6 . In the MS scheme the gluon is more important for jets at 
high x at NNLO than at NLO because the high- a; quarks are smaller. 

Hence, we have included both new data and new theoretical correc- 
tions in our global analysis. The NNLO fit improves on that at NLO. For 
both the value of ag(Mf ) creeps upwards. The NNLO procedure is essen- 
tially complete and we have a preliminary update of parton distributions. 
There are more new data to be included - HERA jets 17 , updated Tevatron 
high-i^T jets, any further new heavy flavour data from HERA and a full 
treatment of NuTcV di-muon data. We also need some further theoretical 
fine-tuning, and will have fully updated NLO and NNLO partons for the 
LHC complete with uncertainties - both experimental and theoretical. 
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